HLA-G Gene Variability Is Associated with Papillary Thyroid Carcinoma Morbidity and the HLA-G Protein Profile.

Human leukocyte antigen (HLA)-G is an immune checkpoint molecule that is highly expressed in papillary thyroid carcinoma (PTC). The HLA-G gene presents several functional polymorphisms distributed across the coding and regulatory regions (5'URR: 5' upstream regulatory region and 3'UTR: 3' untranslated region) and some of them may impact HLA-G expression and human malignancy. To understand the contribution of the HLA-G genetic background in PTC, we studied the HLA-G gene variability in PTC patients in association with tumor morbidity, HLA-G tissue expression, and plasma soluble (sHLA-G) levels. We evaluated 185 PTC patients and 154 healthy controls. Polymorphic sites defining coding, regulatory and extended haplotypes were characterized by sequencing analyses. HLA-G tissue expression and plasma soluble HLA-G levels were evaluated by immunohistochemistry and ELISA, respectively. Compared to the controls, the G0104a(5'URR)G*01:04:04(coding)UTR-03(3'UTR) extended haplotype was underrepresented in the PTC patients, while G0104a(5'URR)G*01:04:01(coding)UTR-03(3'UTR) was less frequent in patients with metastatic and multifocal tumors. Decreased HLA-G tissue expression and undetectable plasma sHLA-G were associated with the G010102a(5'URR)G*01:01:02:01(coding)UTR-02(3'UTR) extended haplotype. We concluded that the HLA-G variability was associated with PTC development and morbidity, as well as the magnitude of the encoded protein expression at local and systemic levels.

BRAF, TERT and HLA-G Status in the Papillary Thyroid Carcinoma: A Clinicopathological Association Study.

As BRAF, TERT, HLA-G, and microRNAs have been individually associated with papillary thyroid carcinoma (PTC), we aimed to evaluate the individual and collaborative role of these markers in PTC in the same patient cohort. HLA-G and BRAF tumor expression was evaluated by immunohistochemistry. Using molecular methods, BRAFV600E and TERT promoter mutations were evaluated in thyroid fine needle aspirates. MicroRNA tumor profiling was investigated using massively parallel sequencing. We observed strong HLA-G (67.96%) while BRAF (62.43%) staining was observed in PTC specimens. BRAF overexpression was associated with poor response to therapy. The BRAFV600E (52.9%) and TERTC228T (13%) mutations were associated with extrathyroidal extension, advanced-age, and advanced-stage cancer. The TERT rs2853669 CC+TC genotypes (38%) were overrepresented in metastatic tumors. Nine modulated microRNAs targeting the BRAF, TERT, and/or HLA-G genes were observed in PTC and involved with cancer-related signaling pathways. The markers were individually associated with PTC features, emphasizing the synergistic effect of BRAFV600E and TERTC228T; however, their collaborative role on PTC outcome was not fully demonstrated. The differentially expressed miRNAs targeting the BRAF and/or HLA-G genes may explain their increased expression in the tumor milieu.

The + 3010/C single nucleotide polymorphism (rs1710) at the HLA-G 3' untranslated region is associated with a short transcript exhibiting a deletion of 92 nucleotides.

The physiological expression of HLA-G is mainly observed in the placenta, playing an essential role in maternal-fetal tolerance. Among the HLA-G mRNA alternative transcripts, the one lacking 92 bases at the HLA-G 3' untranslated region (3'UTR), the 92bDel transcript, is more stable, is associated with increased HLA-G soluble levels, and was observed in individuals presenting a 14 bp insertion (14 bp+) at the 3'UTR. We investigated the presence of the 92bDel transcript in placenta samples, correlating its expression levels with the HLA-G polymorphisms at the 3'UTR. The 14 bp+ allele correlates with the presence of the 92bDel transcript. However, the polymorphism triggering this alternative splicing is the + 3010/C allele (rs1710, allele C). Most 14 bp+ haplotypes (UTR-2/-5/-7) present allele + 3010/C. However, 14 bp- haplotypes such as UTR-3 are also associated with + 3010/C, and the 92bDel transcript can be detected in homozygous samples for the 14 bp- allele carrying at least one copy of UTR-3. The UTR-3 haplotype is associated with alleles G*01:04 and the HLA-G lineage HG0104, which is a high-expressing lineage. The only HLA-G lineage that is not likely to produce this transcript is HG010101, associated with the + 3010/G allele. This functional difference may be advantageous, considering the high worldwide frequency of the HG010101 lineage. Therefore, HLA-G lineages are functionally distinct regarding the 92bDel transcript expression, and the 3010/C allele triggers the alternative splicing that produces this shorter and more stable transcript.

Computational and atomistic studies applied to the understanding of the structural and behavioral features of the immune checkpoint HLA-G molecule and gene.

We took advantage of the increasingly evolving approaches for in silico studies concerning protein structures, protein molecular dynamics (MD), protein-protein and protein-DNA docking to evaluate: (i) the structure and MD characteristics of the HLA-G well-recognized isoforms, (ii) the impact of missense mutations at HLA-G receptor genes (LILRB1/2), and (iii) the differential binding of the hypoxia-inducible factor 1 (HIF1) to hypoxia-responsive elements (HRE) at the HLA-G gene. Besides reviewing these topics, they were revisited including the following novel results: (i) the HLA-G6 isoforms were unstable docked or not with ß2-microglobulin or peptide, (ii) missense mutations at LILRB1/2 genes, exchanging amino acids at the intracellular domain, particularly those located within and around the ITIM motifs, may impact the HLA-G binding strength, and (iii) HREs motifs at the HLA-G promoter or exon 2 regions exhibiting a guanine at their third position present a higher affinity for HIF1 when compared to an adenine at the same position. These data shed some light into the functional aspects of HLA-G, particularly how polymorphisms may influence the role of the molecule. Computational and atomistic studies have provided alternative tools for experimental physical methodologies, which are time-consuming, expensive, demanding large quantities of purified proteins, and exhibit low output.

The Severity of COVID-19 Affects the Plasma Soluble Levels of the Immune Checkpoint HLA-G Molecule.

The non-classical histocompatibility antigen G (HLA-G) is an immune checkpoint molecule that has been implicated in viral disorders. We evaluated the plasma soluble HLA-G (sHLA-G) in 239 individuals, arranged in COVID-19 patients (n = 189) followed up at home or in a hospital, and in healthy controls (n = 50). Increased levels of sHLA-G were observed in COVID-19 patients irrespective of the facility care, gender, age, and the presence of comorbidities. Compared with controls, the sHLA-G levels increased as far as disease severity progressed; however, the levels decreased in critically ill patients, suggesting an immune exhaustion phenomenon. Notably, sHLA-G exhibited a positive correlation with other mediators currently observed in the acute phase of the disease, including IL-6, IL-8 and IL-10. Although sHLA-G levels may be associated with an acute biomarker of COVID-19, the increased levels alone were not associated with disease severity or mortality due to COVID-19. Whether the SARS-CoV-2 per se or the innate/adaptive immune response against the virus is responsible for the increased levels of sHLA-G are questions that need to be further addressed.

The AIRE G228W mutation disturbs the interaction of AIRE with its partner molecule SIRT1.

The autoimmune regulator (AIRE) protein functions as a tetramer, interacting with partner proteins to form the "AIRE complex," which relieves RNA Pol II stalling in the chromatin of medullary thymic epithelial cells (mTECs). AIRE is the primary mTEC transcriptional controller, promoting the expression of a large set of peripheral tissue antigen genes implicated in the negative selection of self-reactive thymocytes. Under normal conditions, the SIRT1 protein temporarily interacts with AIRE and deacetylates K residues of the AIRE SAND domain. Once the AIRE SAND domain is deacetylated, the binding with SIRT1 is undone, allowing the AIRE complex to proceed downstream with the RNA Pol II to the elongation phase of transcription. Considering that the in silico and in vitro binding of the AIRE SAND domain with SIRT1 provides a powerful model system for studying the dominant SAND G228W mutation mechanism, which causes the autoimmune polyglandular syndrome-1, we integrated computational molecular modeling, docking, dynamics between the whole SAND domain with SIRT1, and surface plasmon resonance using a peptide harboring the 211 to 230 residues of the SAND domain, to compare the structure and energetics of binding/release between AIRE G228 (wild-type) and W228 (mutant) SAND domain to SIRT1. We observed that the G228W mutation in the SAND domain negatively influences the AIRE-SIRT1 interaction. The disturbed interaction might cause a disruption in the binding of the AIRE SAND domain with the SIRT1 catalytic site, impairing the AIRE complex to proceed downstream with RNA Pol II.

Placental Malaria is Associated with Higher LILRB2 Expression in Monocyte Subsets and Lower Anti-Malarial IgG Antibodies During Infancy.

Background: Placental malaria (PM) is associated with a higher susceptibility of infants to Plasmodium falciparum (Pf) malaria. A hypothesis of immune tolerance has been suggested but no clear explanation has been provided so far. Our goal was to investigate the involvement of inhibitory receptors LILRB1 and LILRB2, known to drive immune evasion upon ligation with pathogen and/or host ligands, in PM-induced immune tolerance. Method: Infants of women with or without PM were enrolled in Allada, southern Benin, and followed-up for 24 months. Antibodies with specificity for five blood stage parasite antigens were quantified by ELISA, and the frequency of immune cell subsets was quantified by flow cytometry. LILRB1 or LILRB2 expression was assessed on cells collected at 18 and 24 months of age. Findings: Infants born to women with PM had a higher risk of developing symptomatic malaria than those born to women without PM (IRR=1.53, p=0.040), and such infants displayed a lower frequency of non-classical monocytes (OR=0.74, p=0.01) that overexpressed LILRB2 (OR=1.36, p=0.002). Moreover, infants born to women with PM had lower levels of cytophilic IgG and higher levels of IL-10 during active infection. Interpretation: Modulation of IgG and IL-10 levels could impair monocyte functions (opsonisation/phagocytosis) in infants born to women with PM, possibly contributing to their higher susceptibility to malaria. The long-lasting effect of PM on infants' monocytes was notable, raising questions about the capacity of ligands such as Rifins or HLA-I molecules to bind to LILRB1 and LILRB2 and to modulate immune responses, and about the reprogramming of neonatal monocytes/macrophages.

Bullous pemphigoid and milia: prevalence and clinical laboratory findings in a Brazilian sample

Abstract Background: Bullous pemphigoid (BP) associated with milia lesions has been increasingly reported, but its prevalence has not been reported in a Brazilian BP population yet. Objectives: To describe the occurrence and clinical-laboratorial findings of BP-milia association in a southeastern Brazilian sample. Methods: A descriptive study based on the medical charts of 102 BP patients was accomplished. Clinical and laboratory data of BP-milia patients were compiled. Total serum IgE measurements, immunoblot assays based on basement membrane zone antigens, and HLA-DQ alleles typing were performed. Results: Milia was evident in 8 (7.8%) BP patients, five males, aged between 46 and 88 years. Increased total IgE levels were determined in 7 (87.5%) of the eight patients. In five of eight patients, immunoblotting showed IgG reactivity against the BP180-NC16a domain but not against collagen VII or laminin-332; it also revealed reactivity against the BP180 C-terminal domain or LAD-1, or both in four of them. The HLA-DQB1*03:01 and HLA-DQA1*05:05 alleles were identified in three of five BP-milia patients. Moreover, three of five cases presented the HLA-DQB1*06 allelic group. Study limitations: HLA determination was performed in five patients. Conclusions: Milia formation in BP patients seems to be less uncommon than previously admitted. Laboratory data revealed increased IgE; autoantibodies against the BP180 C-terminal domain or LAD-1, or both; and the HLA-DQB1*06 allelic group, described for the BP-milia association. Careful determination of antibodies against basement membrane zone molecules and HLA characterization in different populations may provide further insights into this association. © 2022 Sociedade Brasileira de Dermatologia. Published by Elsevier España, S.L.U. This is an open access article under the CC BY license (https://creativecommons.org/licenses/by/4.0/).

Genetic diversity of the LILRB1 and LILRB2 coding regions in an admixed Brazilian population sample.

Leukocyte immunoglobulin (Ig)-like receptors (LILR) LILRB1 and LILRB2 may play a pivotal role in maintaining self-tolerance and modulating the immune response through interaction with classical and nonclassical HLA molecules. Although both diversity and natural selection patterns over HLA genes have been extensively evaluated, little information is available concerning the genetic diversity and selection signatures on the LILRB1/2 regions. Therefore, we identified the LILRB1/2 genetic diversity using next-generation sequencing in a population sample from São Paulo State, Brazil. We identified 58 LILRB1 Single Nucleotide Variants (SNVs), which gave rise to 13 haplotypes, and 41 LILRB2 SNVs arranged into 11 haplotypes. Although we may not exclude as a possible effect of population structure, we found evidence of either positive or purifying selection on LILRB1/2 coding regions. Some residues in both proteins showed to be under the effect of positive selection, suggesting that amino acid replacements in these proteins resulted in beneficial functional changes. Finally, we have revealed that allelic variation (six and five amino acid exchanges in LILRB1 and LILRB2, respectively) affects the structure and/or stability of both molecules. Nonetheless, LILRB2 has shown higher average stability, with no D1/D2 residue affecting protein structure. Overall, our findings demonstrate that LILRB1 and LILRB2 are as polymorphic as HLA class Ib genes and provide strong evidence supporting the directional selection regime hypothesis.

Matrix Metalloproteinases on Severe COVID-19 Lung Disease Pathogenesis: Cooperative Actions of MMP-8/MMP-2 Axis on Immune Response through HLA-G Shedding and Oxidative Stress.

Patients with COVID-19 predominantly have a respiratory tract infection and acute lung failure is the most severe complication. While the molecular basis of SARS-CoV-2 immunopathology is still unknown, it is well established that lung infection is associated with hyper-inflammation and tissue damage. Matrix metalloproteinases (MMPs) contribute to tissue destruction in many pathological situations, and the activity of MMPs in the lung leads to the release of bioactive mediators with inflammatory properties. We sought to characterize a scenario in which MMPs could influence the lung pathogenesis of COVID-19. Although we observed high diversity of MMPs in lung tissue from COVID-19 patients by proteomics, we specified the expression and enzyme activity of MMP-2 in tracheal-aspirate fluid (TAF) samples from intubated COVID-19 and non-COVID-19 patients. Moreover, the expression of MMP-8 was positively correlated with MMP-2 levels and possible shedding of the immunosuppression mediator sHLA-G and sTREM-1. Together, overexpression of the MMP-2/MMP-8 axis, in addition to neutrophil infiltration and products, such as reactive oxygen species (ROS), increased lipid peroxidation that could promote intensive destruction of lung tissue in severe COVID-19. Thus, the inhibition of MMPs can be a novel target and promising treatment strategy in severe COVID-19.

Bullous pemphigoid and milia: prevalence and clinical laboratory findings in a Brazilian sample.

BACKGROUND: Bullous pemphigoid (BP) associated with milia lesions has been increasingly reported, but its prevalence has not been reported in a Brazilian BP population yet. OBJECTIVES: To describe the occurrence and clinical-laboratorial findings of BP-milia association in a southeastern Brazilian sample. METHODS: A descriptive study based on the medical charts of 102 BP patients was accomplished. Clinical and laboratory data of BP-milia patients were compiled. Total serum IgE measurements, immunoblot assays based on basement membrane zone antigens, and HLA-DQ alleles typing were performed. RESULTS: Milia was evident in 8 (7.8%) BP patients, five males, aged between 46 and 88 years. Increased total IgE levels were determined in 7 (87.5%) of the eight patients. In five of eight patients, immunoblotting showed IgG reactivity against the BP180-NC16a domain but not against collagen VII or laminin-332; it also revealed reactivity against the BP180 C-terminal domain or LAD-1, or both in four of them. The HLA-DQB1*03:01 and HLA-DQA1*05:05 alleles were identified in three of five BP-milia patients. Moreover, three of five cases presented the HLA-DQB1*06 allelic group. STUDY LIMITATIONS: HLA determination was performed in five patients. CONCLUSIONS: Milia formation in BP patients seems to be less uncommon than previously admitted. Laboratory data revealed increased IgE; autoantibodies against the BP180 C-terminal domain or LAD-1, or both; and the HLA-DQB1*06 allelic group, described for the BP-milia association. Careful determination of antibodies against basement membrane zone molecules and HLA characterization in different populations may provide further insights into this association.

The absence of the autoimmune regulator gene (AIRE) impairs the three-dimensional structure of medullary thymic epithelial cell spheroids.

BACKGROUND: Besides controlling the expression of peripheral tissue antigens, the autoimmune regulator (AIRE) gene also regulates the expression of adhesion genes in medullary thymic epithelial cells (mTECs), an essential process for mTEC-thymocyte interaction for triggering the negative selection in the thymus. For these processes to occur, it is necessary that the medulla compartment forms an adequate three-dimensional (3D) architecture, preserving the thymic medulla. Previous studies have shown that AIRE knockout (KO) mice have a small and disorganized thymic medulla; however, whether AIRE influences the mTEC-mTEC interaction in the maintenance of the 3D structure has been little explored. Considering that AIRE controls cell adhesion genes, we hypothesized that this gene affects 3D mTEC-mTEC interaction. To test this, we constructed an in vitro model system for mTEC spheroid formation, in which cells adhere to each other, establishing a 3D structure. RESULTS: The comparisons between AIRE wild type (AIREWT) and AIRE KO (AIRE-/-) 3D mTEC spheroid formation showed that the absence of AIRE: i) disorganizes the 3D structure of mTEC spheroids, ii) increases the proportion of cells at the G0/G1 phase of the cell cycle, iii) increases the rate of mTEC apoptosis, iv) decreases the strength of mTEC-mTEC adhesion, v) promotes a differential regulation of mTEC classical surface markers, and vi) modulates genes encoding adhesion and other molecules. CONCLUSIONS: Overall, the results show that AIRE influences the 3D structuring of mTECs when these cells begin the spheroid formation through controlling cell adhesion genes.

Risk factors for mortality from sepsis in an intensive care unit in Ecuador: A prospective study.

ABSTRACT: To investigate risk factors for mortality from sepsis in an intensive care unit (ICU) in Quito-Ecuador and their association to adherence to Surviving Sepsis Campaign recommendations.Prospective cohort study of patients with severe sepsis/septic shock admitted to the ICU of a public Ecuadorian hospital from March, 2018 to March, 2019. Demographic, clinical, treatment, and outcome data were collected from patients' health records. Patients were divided into 2 groups according to ICU survival or death. Log-binomial regression models were used to identify risk factors for mortality.In total, 154 patients were included. Patients who died in the ICU (n = 42; 27.3%) had higher sequential organ failure assessment score (median 11.5 vs 9; P<.01), more organ dysfunction (median 4 vs 3; P<.0001), and received greater volumes of fluid resuscitation in the first 6 hours (median 800 vs 600 mL; P = .01). Dysfunction of > 2 organs was a risk factor for mortality (relative risks [RR] 3.80, 95% CI 1.33-10.86), while successful early resuscitation (RR 0.32, 95% CI 0.15-0.70), successful empirical antibiotic treatment (RR 0.38, 95%CI 0.18-0.82), and antibiotic de-escalation (RR 0.28, 95%CI 0.13-0.61) were protective factors.Dysfunction of >2 organs was a risk factor for mortality from sepsis while successful early resuscitation and appropriate antibiotic treatment were protective.

Transcript Expression Profiles and MicroRNA Regulation Indicate an Upregulation of Processes Linked to Oxidative Stress, DNA Repair, Cell Death, and Inflammation in Type 1 Diabetes Mellitus Patients.

Type 1 diabetes (T1D) arises from autoimmune-mediated destruction of insulin-producing ß-cells leading to impaired insulin secretion and hyperglycemia. T1D is accompanied by DNA damage, oxidative stress, and inflammation, although there is still scarce information about the oxidative stress response and DNA repair in T1D pathogenesis. We used the microarray method to assess mRNA expression profiles in peripheral blood mononuclear cells (PBMCs) of 19 T1D patients compared to 11 controls and identify mRNA targets of microRNAs that were previously reported for T1D patients. We found 277 differentially expressed genes (220 upregulated and 57 downregulated) in T1D patients compared to controls. Analysis by gene sets (GSA and GSEA) showed an upregulation of processes linked to ROS generation, oxidative stress, inflammation, cell death, ER stress, and DNA repair in T1D patients. Besides, genes related to oxidative stress responses and DNA repair (PTGS2, ATF3, FOSB, DUSP1, and TNFAIP3) were found to be targets of four microRNAs (hsa-miR-101, hsa-miR148a, hsa-miR-27b, and hsa-miR-424). The expression levels of these mRNAs and microRNAs were confirmed by qRT-PCR. Therefore, the present study on differential expression profiles indicates relevant biological functions related to oxidative stress response, DNA repair, inflammation, and apoptosis in PBMCs of T1D patients relative to controls. We also report new insights regarding microRNA-mRNA interactions, which may play important roles in the T1D pathogenesis.

Lenvatinib Plus Anti-PD-1 Combination Therapy for Advanced Cancers: Defining Mechanisms of Resistance in an Inducible Transgenic Model of Thyroid Cancer.

Background: Combination therapy with lenvatinib plus programmed death-1 (PD-1) immune checkpoint blockades (ICBs) is under investigation in many solid tumors, including thyroid cancer. Lenvatinib is known to reduce angiogenesis and may overturn the immunosuppressive effects of vascular endothelial growth factor in the tumor microenvironment. Previous studies investigating the effects of VEGF receptor inhibition on the immune response were performed in rapidly growing tumor models where immune equilibrium is not established before treatment. We hypothesize that physiologically relevant preclinical models are necessary to define mechanisms of resistance to immune-targeted combination therapies. Methods: We utilized the TPO-CreER/BrafV600E/wt/Trp53Δex2-10/Δex2-10 inducible transgenic model of advanced thyroid cancer to investigate lenvatinib treatment in the context of an anti-PD-1 ICB. Following tumor establishment, 3.5 months postinduction, mice were treated with high- (10 mg/kg) or low-dose (2 mg/kg) lenvatinib, anti-PD-1, or combination of lenvatinib with anti-PD-1. Tumor volume and lung metastases were assessed in each group. Immune infiltrate was characterized by flow cytometry and immunohistochemistry, and TCRß sequencing was performed to further investigate the T cell response. Results: Both low- and high-dose lenvatinib reduced tumor volume, while anti-PD-1 had no effect, alone or in combination. Although both low- and high-dose lenvatinib reduced vascular density, low-dose lenvatinib was superior in controlling tumor size. Lung metastases and survival were not improved with therapy despite the effects of lenvatinib on primary tumor size. Low-dose lenvatinib treatment led to a subtle reduction in the dominant Ly6G+CD11b+ myeloid cell population and was associated with increased CD4+ T cell infiltrate and enrichment in 4-1BB+ and granzyme B+ CD4+ T cells and FoxP3+ regulatory T cells. Polyclonal T cell expansion was evident in the majority of mice, suggesting that a tumor-specific T cell response was generated. Conclusions: The effects of lenvatinib on the immune response were most pronounced in mice treated with low-dose lenvatinib, suggesting that dose should be considered in clinical application. While the immune-modulating potential of lenvatinib is encouraging, alterations in the immune milieu and T cell activation status were insufficient to sustain durable tumor regression, even with added anti-PD-1. Additional studies are necessary to develop more effective combination approaches in low-mutation burden tumors, such as thyroid cancer.

Structural Characterization of the Interaction of Hypoxia Inducible Factor-1 with Its Hypoxia Responsive Element at the -964G > A Variation Site of the HLA-G Promoter Region.

Human Antigen Leukocyte-G (HLA-G) gene encodes an immune checkpoint molecule that has restricted tissue expression in physiological conditions; however, the gene may be induced in hypoxic conditions by the interaction with the hypoxia inducible factor-1 (HIF1). Hypoxia regulatory elements (HRE) located at the HLA-G promoter region and at exon 2 are the major HIF1 target sites. Since the G allele of the -964G > A transversion induces higher HLA-G expression when compared to the A allele in hypoxic conditions, here we analyzed HIF1-HRE complex interaction at the pair-atom level considering both -964G > A polymorphism alleles. Mouse HIF2 dimer crystal (Protein Data Bank ID: 4ZPK) was used as template to perform homology modelling of human HIF1 quaternary structure using MODELLER v9.14. Two 3D DNA structures were built from 5'GCRTG'3 HRE sequence containing the -964G/A alleles using x3DNA. Protein-DNA docking was performed using the HADDOCK v2.4 server, and non-covalent bonds were computed by DNAproDB server. Molecular dynamic simulation was carried out per 200 ns, using Gromacs v.2019. HIF1 binding in the HRE containing -964G allele results in more hydrogen bonds and van der Waals contact formation than HRE with -964A allele. Protein-DNA complex trajectory analysis revealed that HIF1-HRE-964G complex is more stable. In conclusion, HIF1 binds in a more stable and specific manner at the HRE with G allele.

HLA-G genetic diversity and evolutive aspects in worldwide populations.

HLA-G is a promiscuous immune checkpoint molecule. The HLA-G gene presents substantial nucleotide variability in its regulatory regions. However, it encodes a limited number of proteins compared to classical HLA class I genes. We characterized the HLA-G genetic variability in 4640 individuals from 88 different population samples across the globe by using a state-of-the-art method to characterize polymorphisms and haplotypes from high-coverage next-generation sequencing data. We also provide insights regarding the HLA-G genetic diversity and a resource for future studies evaluating HLA-G polymorphisms in different populations and association studies. Despite the great haplotype variability, we demonstrated that: (1) most of the HLA-G polymorphisms are in introns and regulatory sequences, and these are the sites with evidence of balancing selection, (2) linkage disequilibrium is high throughout the gene, extending up to HLA-A, (3) there are few proteins frequently observed in worldwide populations, with lack of variation in residues associated with major HLA-G biological properties (dimer formation, interaction with leukocyte receptors). These observations corroborate the role of HLA-G as an immune checkpoint molecule rather than as an antigen-presenting molecule. Understanding HLA-G variability across populations is relevant for disease association and functional studies.

Peripheral spectrum neurological disorder after arbovirus infection is associated with HLA-F variants among Northeastern Brazilians.

INTRODUCTION: Non-classical class I human leukocyte antigens (HLA) molecules are known to modulate the function of cytotoxic cells (NK and T CD8+) during viral infection by interacting with inhibitory/activating receptors. However, little is known about the HLA-E/-F genetic variability on arbovirus infections. METHODS: We evaluated by massive parallel sequencing the full HLA-E/-F genetic diversity among patients infected during the arbovirus (ZIKV, DENV, and CHIKV) outbreak leading to a broad range of neurological complications in the Brazilian State of Pernambuco. In parallel, healthy blood donors from the same area were also studied. Plink and R software were used for genetic association study. To limit the false-positive results and enhance the reliability of the results, we adopted P-values <0.01 as significant levels. RESULTS: Compared to controls, the HLA-F alleles: -1610 C (rs17875375), +1383 G (rs17178385), and +3537 A (rs17875384), all in complete linkage disequilibrium with each other (r2 = 1), were overrepresented in patients presenting peripheral spectrum disorders (PSD). The HLA-F*Distal-D haplotype that harbored the -1610 C allele exhibited a trend increase in PSD group. No associations were found for HLA-E. CONCLUSIONS: Our findings showed that the HLA-F genetic background seems to be more important than HLA-E on the susceptibility to PSD complications.

Human leukocyte antigen (HLA)-F and -G gene polymorphisms and haplotypes are associated with malaria susceptibility in the Beninese Toffin children.

BACKGROUND: Little attention has been devoted to the role of the immunoregulatory HLA-E/-F/-G genes in malaria. We evaluated the entire HLA-E/-F/-G variability in Beninese children highly exposed to Plasmodium falciparum (P.f.) malaria. METHODS: 154 unrelated children were followed-up for six months and evaluated for the presence and number of malaria episodes. HLA-E/-F/-G genes were genotyped using massively parallel sequencing. Anti P.f. antibodies were evaluated using ELISA. RESULTS: Children carrying the G allele at HLA-F (-1499,rs183540921) showed increased P.f. asymptomatic/symptomatic ratio, suggesting that these children experienced more asymptomatic P.f. episodes than symptomatic one. Children carrying HLA-G-UTR-03 haplotype exhibited increased risk for symptomatic P.f. episodes and showed lower IgG2 response against P.f. GLURP-R2 when compared to the non-carriers. No associations were observed for the HLA-E gene. CONCLUSION: HLA-F associations may be related to the differential expression profiles of the encoded immunomodulatory molecules, and the regulatory sites at the HLA-G 3'UTR may be associated to posttranscriptional regulation of HLA-G and to host humoral response against P.f.